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V0332+53 in the outburst of 2004-2005: luminosity 
dependence of the cyclotron line and pulse profile 
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ABSTRACT 

We present results of observations of the transient X-ray pulsar V0332-I-53 performed 
during a very powerful outburst in Dec, 2004 - Feb, 2005 with the INTEGRAL and 
RXTE observatories in a wide (3 — 100 keV) energy band. A cyclotron resonance 
scattering line at an energy of ~ 26 keV has been detected in the source spectrum 
together with its two higher harmonics at ~ 50 and ~ 73 keV, respectively. We show 
that the energy of the line is not constant but linearly changes with the source lumi- 
nosity. Strong pulse profile variations, especially near the cyclotron line, are revealed 
for different levels of the source intensity. We discuss the obtained results in terms of 
the theoretical models of X-ray pulsars. 
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1 INTRODUCTION 

, The transient X-ray pulsar V0 332-I-53 was discovered by 
■ the Vela 5B observatory in 1973 jXerrel fc Predhorskvll97a) 
during an outburst when its intensity reached ~ 1.4 Crab 
r in the 3 — 12 keV energy band. The outburst lasted about 
three months before the source became undetectable again. 
During observations in Nov 1983 - Jan 1984 with the 
' EXOSAT observatory the pulsar's parameters and orbital 
ones were determined: the pulse period ~ 4.375 s, the or- 
bital period - 34.25 days, the eccentricity - 0.31, and the 
projected somimajor axis of the neutron star - a^^sini ~ 48 
It-s (.Stella ct al. 1985). These authors also mentioned that 
as the source intensity decreased the pulse profile changed 
from double to single peaked, which was accompanied by 
significant hardening of the source spectrum. Later when 
the source was observed by the GINGA observatory the 
cyclotron resonance scattering feature with an energy of 
Ecyc ~ 28.5 ± 0.5 keV was detected in its spectrum. This 
energy corresponds to a magnetic field on the neutron star 
surface of ~ 3 x 1 0^^ G ( Makishima ct al. 1990). Later 
iMihara et al.l (Il998ft reported measurements of two diflfer- 
ent values of the resonance energy for different levels of the 
source intensity. 

The next powerful outburst of the source began at the 
end of 2004 ( Swank ct al. 2004). This outburst had been 
predicted based on the increasing of the optical brightness 
of the norma l companion which reached i ts maximum on 
31 Jan 2004 iCoranskii fc Barsukovall20oil . A preliminary 
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analysis of RXTE observations performed on 24-26 Dec 
2004 showed that beside the absorption feature at an en- 
ergy of 26.34 ± 0.03 keV, there are two additional similar 
features in the source spectrum at energies of 49.1 ± 0.2 
and 74 ± 2 keV. These features were interpreted as the sec- 
ond and third harm onics of the main cyclotron frequency 
(ICoburn et alj|2005ll . The similar results were inferred from 
analysis of the first ^ 100 ks data of the INTEGRAL obser- 
vatory dKrevkenbohm et alj l2005'l. 



INegueruela et aP (^9^ discussed the results of the op- 
tical observations of BQ Cam - the normal companion of the 
X-ray pulsar V0332-1-53. They determined the star's spectral 
class as 08-9Ve and estimated the distance to the system at 
~ 7 kpc. 

In this paper we present results of an analysis of obser- 
vations of V0332-f 53 performed with the INTEGRAL and 
RXTE observatories in Jan-Feb 2005. Our main aim is to 
study variations of the source spectrum and pulse profile 
depending on its intensity. Results of the spectral and tim- 
ing analysis based on RXTE data and in dependent analysis 
of the INT EGRAL data can be found in IPottschmidt et alJ 
(|2005^ and lMowlavi et al.l J2006l) . respectively. The observa- 
tions, instruments and data analysis are described in section 
2. In section 3 we demonstrate for the first time that the 
energy of the cyclotron resonance scattering feature is not 
constant but linearly increases with decreasing luminosity 
of the source. Section 4 is dedicated to a detailed study of 
the pulse profile variations especially near the cyclotron line. 
The obtained results are discussed in section 5. 
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2 OBSERVATIONS AND DATA ANALYSIS 

The X-ray pulsar V0332-I-53 was observed with the IN- 
TEGRAL and RXTE observatories several times at the 
end of 2004 - beginning of 2005. About 130 pointed ob- 
serva tions were carried o ut by the INTEGRAL observa- 
tory dWinkler et al.ll2003D during TOO observations of the 
source with a total exposure ~ 400 ks (these data are 
public a vailable). Data of th e IBIS telescope (the ISGRI 
detector , iLebrun et al.l i2003ft ') and X-ray monitor JEM-X 
jLund et al J 1200311 were used in the analysis. The effective 
energy bands of these instruments are 20 — 200 and 3 — 30 
keV, respectively, which makes it possible to study bright 
astrophysical objects in a wide energy band. 

The image reconstructi on was done with the m ethod 
and software described in iRevnivtsev et all (|20o3). The 
source spectra in a hard energy band (> 20 keV) obtained by 
the IBIS telescope were calculated through reconstruction 
of a large number of images of the source in narrow energy 
channels and flux extraction of the studied source. The re- 
sponse matrix for our research have been done based on the 
standard INTEGRAL matrix (taken from the OSA pack- 
age), but the new arf-file was calculated using the calibra- 
tion observations of the Crab Nebulae. Standard calibration 
tables and procedures (similar to those implemented in OSA 
5.0) were used to correct the energy of ISGRI events having 
different rise-time. The analysis of a large number of cali- 
bration observations for the Crab nebula revealed that this 
method has a systematic error in measuring source fluxes 
about ^ 3%. We included this systematic uncertainty when 
analyzing spectra with the XSPEC package. 

For the timing analysis on time scales of the order of the 
source pulse period and for the analysis of the JEM-X data 
we used the standard OSA 5.0 software package provided by 
the INTEGRAL Science Data Center (http://isdc.unige.ch). 
For the pulse profiles reconstruction at high energies the first 
three steps of OSA (COR, GTI, DEAD) were executed. Af- 
ter this we collected photons with specified energies from 
the detector pixels opened for the source using the tool 
evts-extract. At the next stage, the arrival times of collected 
photons were corrected for the neu tron star orbital m otion 
using the known binary parameters iStella et al Eiil. Such 
a relatively simple approach for collecting photons is justi- 
fiable in our case as the source V0332-I-53 was a single ob- 
ject in the telescope field of view. Thus the detector collects 
only photons from the source and background photons. The 
background can be calculated from those pixels which are 
fully closed for the source by the opaque mask elements. To 
reconstruct the pulse profile at low energies (JEM-X data) 
photons from the whole detector were collected (the tool 
evts_pick). The subsequent timing analysis was done with 
the FTOOLS package. 

Due to a nicety of the considering effects (in particular 
a displacement of the cyclotron line center is about ~ 3 keV 
in the INTEGRAL data) and the using of different types of 
the software on difi'erent stages of the analysis we carried 
out an additional investigation of an accuracy and stability 
of the energy determination for both type of the software. 
The tungsten Kc, line 59.32 keV was choose as a calibration 
line due to its brightness (strictly speaking, in this spectral 
region there is a blend of tungten lines at energies 57.98, 
59.32 etc. keV, but line at 59.32 keV is the most strong). The 
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Figure 1. Dependence of the energy of the tungsten line 59.3 keV, 
measuring by the ISGRI detector, from the observation date for 
special (a) and standard OSA-5.0 (b) softwares. 



detector spectrum was built for both softwares (for OSA 5.0 
it was done after the DEAD level), then the energy of the line 
center was determined and compared with ones from other 
observations and its theoretical value. It was shown that the 
line center position is stable and in a good agreement with 
the theoretical one in the spectral analysis (Fig^). The 
conservative estimation of the systematic uncertainty of the 
energy determination gives ~ 0.1 keV. 

For the standard software OSA 5.0 the line center posi- 
tion is changed from one observation to another and displace 
from the theoretical value (Fig^). Therefore in the follow- 
ing analysis we took it into account and made a correction 
for this effect. 

Apart from INTEGRAL dat a we also use d 
data of the RXTE observatory jBradt et alJ Il993l) : 
simultaneous data of the all-sky monitor ASM 
(http:/ /xte.mit.edu/ASMJc.htmlf), and the spectrome- 
ters PCA and HEXTE (Obs.ID 90427), operating in the 
energy bands 1.3 — 12.2, 3 — 20 and 15 — 250 keV, respec- 
tively. For the RXTE data reduction we used standard 
programs of FTOOLS/LHEASOFT 5.3 package. As for 
the ISGRI detector, we investigated the accuracy of the 
energy determination using the HEXTE spectrometer. In 
this analysis we used the ^^^Am calibration line at 59.6 
keV. To estimate the systematic error of energy measuring 
we determined a typical scatter of measurements of the line 
center, ~ 0.1 keV. This value was added as a systematic 
error to the uncertainties obtained from the spectral 
analysis. 

A journal of observations of the X-ray pulsar V0332-f 53 
is presented in TabllTI The date of observations, correspond- 
ing exposure, observed source flux and luminosity are given. 
The luminosity was calculated assuming a source distance 
of 7 kpc. The obtained luminosity values may be regarded 
as close to the bolometric ones on the assumption that the 
bulk of the energy is released in X-rays. 
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Table 1. Observations of the X-ray pulsar V0332+53 with the INTEGRAL and RXTE observatories in 2004-2005 



Date, MJD 


Exposure, 


Flux^-, 


Luminosity'-', 


(revolution) 


ks xlO" 


" erg s cm 


10''^ erg s~^ 


INTEGRAL data 


53376.5 


(272 rev) 


28.7 


58.3 ±0.9 


34.1 ± 0.5 


53379.2 


(273 rev) 


56.7 


45.6 ± 0.4 


26.7 ±0.2 


53380.3 


(274 rev) 


21.1 


40.7 ±2.4 


23.8 ± 1.4 


53394.0 


(278 rev) 


72.3 


24.9 ± 1.0 


14.5 ± 0.6 


53410.9 


(284 rev) 


149.1 


12.5 ± 0.2 


7.3 ±0.1 


53413.1 


(285 rev) 


15.5 


11.4 ±0.3 


6.7±0.2 


53416.5 


(286 rev) 


15.4 


8.4 ±0.3 


4.9 ±0.2 


53419.1 


(287 rev) 


17.0 


7.0 ±0.3 


3.4 ±0.2 


53422.1 


(288 rev) 


20.1 


3.5 ±0.2 


2.0 ±0.1 


RXTE data 


53367.2 


(90427-01-01-OOG) 


1.5 


82.4 ±2.0 


48.4 ± 1.2 


53368.3 


(90427-01-01-01) 


2.2 


81.6 ± 2.4 


48.0 ± 1.4 


53368.9 


(90427-01-01-02) 


1.8 


79.2 ± 0.5 


46.7 ±0.3 


53369.6 


(90427-01-01-03) 


1.9 


72.9 ± 0.4 


42.9 ± 0.2 


53374.0 


(90427-01-02-00) 


0.9 


68.5 ± 2.1 


40.3 ± 1.2 


53375.0 


(90427-01-02-01) 


0.8 


60.9 ± 2.1 


35.8 ± 1.2 


53376.3 


(90427-01-02-02) 


0.7 


66.0 ±2.6 


38.8 ± 1.5 


53376.7 


(90427-01-02-03) 


2.7 


61.9 ± 1.3 


36.4 ±0.8 


53385.1 


(90427-01-03-01) 


9.9 


43.7 ± 1.5 


25.7 ±0.9 


53385.5 


(90427-01-03-02) 


14.5 


41.5 ± 0.2 


24.4 ±0.1 


53387.0 


(90427-01-03-05) 


12.8 


38.0 ±0.1 


22.3 ± 0.1 


53387.4 


(90427-01-03-06) 


12.1 


36.5 ± 0.1 


21.5 ± 0.1 


53387.9 


(90427-01-03-07) 


9.6 


37.2 ± 0.1 


21.9 ± 0.1 


53388.4 


(90427-01-03-09) 


11.9 


34.5 ± 0.8 


20.3 ± 0.5 


53389.0 


(90427-01-03-11) 


9.9 


34.5 ± 0.4 


20.3 ± 0.2 


53389.3 


(90427-01-03-12) 


9.7 


34.1 ± 1.0 


20.1 ± 0.6 


53413.1 


(90427-01-04-00) 


6.1 


14.7 ±0.6 


8.6 ±0.3 


53413.8 


(90427-01-04-04) 


6.7 


12.6 ± 0.4 


7.4 ±0.2 


53414.3 


(90427-01-04-02) 


11.8 


11.1 ± 0.2 


6.5 ±0.1 


53414.5 


(90427-01-04-03) 


6.8 


11.0 ±0.3 


6.5 ±0.2 


53414.8 


(90427-01-04-05) 


2.6 


10.9 ± 0.2 


6.4 ±0.1 


53416.6 


(90427-01-04-01) 


5.7 


8.9 ±0.2 


5.2 ±0.1 


a 


in the 3-100 keV 


energy band 






b 


in the 3-100 keV energy band, assuming source distance of 7 kpc 



3 LIGHT CURVES 

The daily-averaged light curve of the X-ray pulsar 
V0332+53 obtained by ASM/ RXTE is presented in FigUa,. 
The source intensity was increasing practically linearly for 
~ 30 days, then the source stayed for about 10 days near 
the maximum of its intensity, and after that the measured 
flux decreased to the pre-outburst value in ~ 50 days. The 
last segment of the light curve can be described by an ex- 
ponential decay with a characteristic time r ~ 17 days. A 
compariso n of the 2004-2005 outburst with previous type II 
outbursts iStella et al.l ll9861 shows that it is a typical out- 
burst of this class for V0332+53 in terms of duration and 
maximum of X-ray flux. 

INTEGRAL observations started only several days after 
the outburst maximum due to limitations in the satellite 
orientation relative to the Sun. In Fig|5|D the pulsar light 
curve obtained by the ISGRI detector in the hard (18 — 60 
keV) energy band is shown. Each point corresponds to a 
separate observation with exposure about 2 ks. The flgure 



demonstrates that the source flux in hard X-rays decreased 
from ~ 900 to ^ 100 mCrab during a month and a half. Note 
that the shape of the light curve in hard X-rays is slightly 
different from the one in the soft 1.3—12.2 keV energy band. 



From the analysis of the V0332+53 light curve in the 
soft energy band it was found that the source intensity 
demonstrates variability with an amplitude of about 20% 
(near the outburst maximum). The variability amplitude 
decreased when the source flux decreased and practically 
disappeared by the end of the outburst. The variability with 
the same amplitude was observed with the ISGRI detector 
in the hard energy band too. The 18-60 keV source light 
curve with time resolution of 300 s is shown in Fig|2J; for 
the bright state of the source, when the variability is vis- 
ible most. In the right part of the Figure a local outburst 
with maximum at MJD 53380.3 and total duration about 
4.5 hours is seen. The solid line represents the best-flt ap- 
proximation by a gaussian. 
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Figure 2. Light curves of tlie X-ray pulsar V0332+53 obtained 
during tlie outburst: (a) with the ASM/ RXTE monitor; (b) with 
the IBIS/ INTEGRAL telescope in a hard energy band; (c) with 
the IBIS/ INTEGRAL telescope with a time resolution of 300 s for 
the bright state. Solid line represents the best-fit approximation of 
the small outburst by a gaussian with maximum at MJD 53380.3 
and total width ~ 4.5 h. 



4 SPECTRAL ANALYSIS 

The V0332-I-53 spectrum deserves special attention. This is 
only the second X-ray accreting pulsar (after 4U0115-I-63) 
whose spectrum exhibit not only a cyclotron reso- 
nance sc attering featu r e bu t also its two highe r har- 
monics JCoburn et all 120051: iKrevkenbohm et alj 120051: 
|Pottschniid^^Lf]2005ll . As a whole, the spectrum of the 
pulsar can be well described by a power law with an ex- 
ponential cutoff at high energies, typical for this class of 
objects. 

The significant number of V0332-I-53 observations car- 
ried out with the INTEGRAL and RXTE observatories al- 
lowed us to reconstruct the source spectrum at different 
phases of the outburst and to trace the evolution of its 
parameters. We approximated the obtained spectra with 
a powerlaw model with an exponential cutoff {cutoffpl in 
XSPEC), modified by three absorption lines in the form 



exp 



ycl{E/ Ecycl 



I "cycl 



{E-E,y^iY + al 



where E^yd, Ocyd 
and depth of the line, 



cycl 

and Tcyci ar e the 
respectively 



center, width 
Il99cll . 



Mihara et alJ 

This model describes the source spectrum similarly well as 
the usual model of a powerlaw with a high energy cutoff 
{powerlaw x highecut in XSPEC, White et al. 1983) but 
has fewer parameters. The best-fit value of the cutoff energy 
EciLt for the powerlaw x highecut model is about (5 — 6) 
keV. Given the limited energy range of the used instru- 
ments (down to 3 keV for PGA and 4.5 keV for JEM-X), 



0.01 ^ 




INTEGRAL/JEM-X 
INTEGRAL/IBIS/ISGRI 



5 10 20 40 60 100 

Energy, keV 

Figure 3. Energy spectra of V0332+53 measured with the IN- 
TEGRAL observatory for two states (272 and 284 revolutions). 



it is not possible to put tight constraints on the photon in- 
dex. An iron emission line was added to the model when 
fitting the RXTE data. This feature is not detected by the 
JEM-X monitor partly due to its lower sensitivity compared 
to PGA, and also due to the current uncertainty in the 
JE M-X response matr ix at these energies (see comments in 
Fili ppova et alJ l|200^ for details). 

A typical source spectrum obtained with the INTE- 
GRAL observatory for the bright (272 rev) state is shown 
in Fig|3] The best-fit parameters are summarized below: 



Model parameter Value 



Photon index 
Ec^t, keV 

'^cycl , 1 

Ecycl, 1, keV 
(^cyci,i, keV 

Tcycl,2 

Ecycl, 2, keV 
f^cyci,2, keV 

'^cycl,3 

Ecycl, 3, keV 
cycyci,3, keV 
(d.o.f) 



-0.120 ±0.008 
9.21 ±0.04 
1.91 ±0.02 

25 09+^ °^ 



5.44^ 



1+0.08 
-0.06 

2.12 ±0.03 
49 44+'^-°'' 

».03_Q 23 

1.26 ±0.10 

79 1+0.5 

1.25(136) 



For comparison, the spectrum obtained in the low state 
(284 rev) is shown in the same Figure. 

The cyclotron resonance scattering feature and its sec- 
ond harmonics are clearly visible in both spectra. Despite 
the rapid decrease of the source intensity and its weakness 
at high energies (> 65 keV), inclusion into the model of 
the third harmonics with an energy of ~ 75 — 80 keV leads 
to a significant improvement of the fit (Ax^ = 18 for 3 
d.o.f.). This harmonics is also detected in several next ob- 
servations, but its parameters (width and depth) are rea- 
sonably bounded only in the bright state (till 284 rev) . Fix- 
ing them at the values obtained for the bright state, does 
not allow one to improve the quality of the fit for the low- 
state spectra. Moreover, the determined line widths for these 
observations are too large and strongly affect the determi- 
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Figure 4. The cyclotron line energy dependence on the source lu- 
minosity (3-100 keV). Triangles are INTEGRAL results, squares 
are RXTE ones. 



nation of the parameters of the second harmonics, mak- 
ing this task model dependent. It is necessary to note that 
the fundamental line energy independs on the including to 
the model the third harmonics. It is nessecary to note that 
IPot tschmidt ot al. (2005) used a Gaussian profile for the de- 
scribing of the cyclotron line, therefore obtained by them 
values of the cyclotron energy are slightly differ from our 
ones. The same situation was discussed bv lNakaiima et al] 
l|200fih for 4U0115-^634. 

Our analysis showed that the model described approx- 
imates the source spectrum well during the whole outburst 
both for INTEGRAL and RXTE data. The behaviour of 
the cyclotron line is of greatest interest because it is confi- 
dently detected during the entire outburst and its parame- 
ters are model independent. The line energy dependence on 
the source luminosity obtained from INTEGRAL and RXTE 
data is shown in Fig|l]by dark triangles and squares, respec- 
tively. The uncertainties of the HEXTE results are slightly 
higher than those for ISGRI as the HEXTE exposures are 
shorter (see Table 0. The measurements of both observato- 
ries are in good mutual agreement and fall on a straight line, 
i.e. the cyclotron line energy increases linearly with decreas- 
ing source luminosity. The formal fitting of this dependence 
with a linear relation gives Ecyci,i — — O.IOL37 -I- 28.97 keV, 
where L37 - luminosity in units of 10''^ erg s~^. Believing 
that for low luminosities the emission come practically from 
the neutron star surface (see section 6.1) we can estimate 
the magnetic field on the surface 



Bns 



28.97 



2GM_NS 



11.6 



where -Rjvs and Mms ~ are the neutron star radius and 
mass, respectively. 

As was mentioned above, the energy of the second har- 
monics is not reasonably determined for all observations 
and depends on the inclusion of the third harmonics in the 
model. To avoid possible contamination by this component 



Figure 5. The same as in Fig|l] but for the second harmon- 
ics. Open triangles represents INTEGRAL measurements for the 
broad energy band (till 110 keV) with inclusion of the third har- 
monics (sec text). 



in determining of the parameters of the second harmonics, 
we restricted the considered energy band by 65 keV and fit- 
ted the spectra by the same model with two absorption lines. 
The variation of the second harmonics with source luminos- 
ity is shown in Fig|^by dark triangles and squares similarly 
to FigH It can be seen that although the scatter is slightly 
larger than for the first harmonics, the overall tendency of 
the line energy increasing with decreasing luminosity holds 
in this case too. Formal fitting of this dependence with a 
linear relation gives Ecyci,2 oc — O.O8L37. The energy of the 
second harmonics was also determined by analyzing the IN- 
TEGRAL data in the broad energy band (up to 110 keV) 
when the third harmonics was added to the model (open 
triangles in Fig|^. The energies obtained in both analyses 
differ from each other, especially for the observations with 
low luminosities. Moreover, the ISGRI measurements lie be- 
low the near simultaneous HEXTE ones. This fact is most 
likely connected with the larger exposures of the ISGRI de- 
tector - its data at high energies provide better statistics and 
the third harmonics affects the determination of the param- 
eters of the second harmonics even in the case of the energy 
band truncated at 65 keV. If we exclude from the consid- 
eration the INTEGRAL results, then formal fitting of the 
HEXTE measurements results in Ecyci,2 oc — Q.IL37, which 
is the same as was obtained earlier for the main harmonics. 



5 PULSE PROFILE 

Because of the high intensity of the pulsar emission we suc- 
ceeded in studying its pulse profile dependence on the energy 
band and luminosity. As the background doesn't affect on 
the pulse profile shape, which is studied in the paper, it was 
not subtracted in the following analysis. The most charac- 
teristic source pulse profiles obtained with the INTEGRAL 
observatory in different energy bands are shown in FigElfor 
two observations (272 and 284 revolutions). The observed 
peculiarities in the pulse profile behavior can be divided in 
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V 0332+53 INTEGRAL/ISGRI + .lEM-X 



INTKCRM./JKM-X [reva7a) ^ V033a + 53 INTEGRAL/JEM -X (rev.264) V0^3a + 53 





Figure 6. V0332+53 pulse profiles obtained with tlie JEM-X 
monitor and IBIS telescope of the INTEGRAL observatory in 
different energy bands for different source intensities. The pulse 
phase zero for each observation are chosen independently. 



two main groups: an asymmetrical evolution of the double- 
peaked profile in wide energy bands and its drastic changes 
near the main harmonics of the cyclotron line. Below we 
describe these effects in detail. 

In the brightest state (272 rev) the pulse profile has a 
sinusoidally-like double-peaked shape with little prevalence 
of the second peak at low energies (3 — 6 keV). As the pho- 
tons energy grows the relative contribution of the first peak 
increases and exceeds the second peak intensity in the 10— 15 
keV energy band. The intensity of the first peak continues to 
grow when increasing the energy and becomes significantly 
larger than the second one in the 30 — 50 keV channel. No 
significant movements of the peaks depending on the energy 
band occur. 

As the pulsar intensity diminishes, significant changes of 
the pulse profile occur in soft energy bands (JEM-X data): 
in the 3 — 6 keV band it becomes nearly single-peaked in 
the 274 revolutions; the relative intensity of the first peak 
is very small in 3 — 6 and 6 — 10 keV bands in 278 and 
284 revolutions. In high-energy channels (ISGRI data) the 
pulse profile evolves similarly as was described above for 
the 272 revolution. In the following observations (beginning 
from the 284 revolution when the source luminosity dropped 
to - 7.3x10^'' ergs"^) the main described tendencies persist 
at soft energies, but at energies of order and above the cy- 
clotron one a drastic change in the pulse profile takes place: 
the intra-peak gap disappears and the double-peaked profile 
transforms into an asymmetrical single-peaked one. When 



the energy increases, the phase of the main minimum is dis- 
placed on ~ 0.5 (Figl^. The described features remain ap- 
proximately the same in wide energy channels for the sub- 
sequent observations, i.e. at lower luminosities. 

A detailed study of pulse profiles was carried out in 
narrow energy channels in the most interesting region - 
around the main harmonics of the cyclotron line (FigEJ. 
Such a study was performed for four INTEGRAL observa- 
tions (revolutions 272, 278, 284 and 286) and for four RXTE 
observations which are close to INTEGRAL ones (90427-01- 
02-03, 90427-01-03- 11+12, 90427-01-04-00 and 90427-01-04, 
respectively). The last was done to obtain an independent 
confirmation of the INTEGRAL results. In Fig|7]pulse pro- 
files for these INTEGRAL observations are drawn. The en- 
ergy channels were chosen to divide the cyclotron line in 
half, their width was approximately the same and about a 
half of the cyclotron line width (i.e. two central channels cor- 
respond to the lower and upper wings of the line). Here we 
took into account that the energy of the line center changes 
with the source luminosity. 

For the brightest state (272 rev), the passage through 
the line center does not affect the pulse profile and all the 
tendencies described above for wide channels remain. In the 
subsequent observations (273-278 revolutions) some changes 
occur in the relative intensity of the peaks and their po- 
sitions, but the pulse profile remains double-peaked. With 
a further decrease of the luminosity to — 7.3 x 10"^^ erg 
s~^ (284 revolution) the profile became asymmetrical single- 
peaked at energies below the cyclotron line with a drastic 
transition to the double-peaked shape above the line energy. 
Unlike in the 278 revolution the displacement of the profile 
by half a period, expressed in the displacement of the main 
minimum, occurs not in the upper wing of the line but in 
the next energy channel (Fig|7J. A similar picture remains 
when the pulsar luminosity is decreased to ~ 4.9 x 10"^^ 
erg s~^ (286 rev). In the next revolution (the luminosity is 
— 3.4 X 10''^ erg s"'^) the pulse profile at energies below 
the line returned to the double-peaked shape, although the 
single-peaked one holds in the lower wing. The displacement 
of the profile by half a period is retained in the last channel 
too. It is necessary to note that due to the source intensity 
decreasing the statistics in the last observations is not suffi- 
cient for detailed analysis of the pulse profile structure and 
we can only consider their common characteristics. 

The results of a similar analysis performed on the 
RXTE /REXTE data (bottom panels on FigQ fully con- 
firm the conclusions about the behavior of the pulse profile 
drawn above based on the INTEGRAL/lSGRl data. 



6 DISCUSSION 

The transient X-ray pulsar V0332-I-53 demonstrates pow- 
erful outbursts in which its intensity exceeds 1 Crab. The 
source is a member of a high mass X-ray binary sys- 
tem with a companion star (BQ Cam) that belongs to 
the class of Be st ars. According to current ideas (e.g. 
[Okazaki fc Negueruela (2001 ) ) such objects represent quick 
rotating stars with a dense, but radially slow, stellar wind. 
This wind forms the so-called equatorial disc around the 
star, the size and the presence of the disc are not perma- 
nent. Most X-ray sources in binary systems with Be stars 
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Figure 7. Evolution of the source pulse profile near the cyclotron line with the source intensity. YBIS/ INTEGRAL data (upper panels), 
HEXTE/RXTE data (bottom panels). The pulse phase zero for each observation are chosen independently. 



are transients demonstrating outbursting activity. This is 
presumably connected with the evolution of the normal star. 
Matter from the equatorial disc is captured and an accretion 
disc around the relativistic object is formed. 

The system V0332+53/BQ Cam is an obvio us example 
of the picture described above. As was shown bv lGoranskiil 
i200ll) . a significant brightening of the normal star in the op- 
tical waveband preceded previous X-ray outbursts. It is be- 
lieved that this brightening is associated with the formation 
and s ubsequent ejectio n of the ambient matter (equatorial 
disc) iGoranskiil 1200 iT) . The outburst that started in Dec, 
2004 was no exception and had been predicted based on an 
increase in the brightness of the optica l star at the beginning 
of 2004 iGoranskii fc B arsukova'2004^ . Such a large (several 
hundred days) delay between the op tical and X-ray light is 
typical for this source lCoranskiil i200ll) and is possibly assco- 
ciated with the necessity to accumulate in the accretion disc 
the sufficient amount of matter to begin accretion. The idea 
that accretion during outbursts proceeds through the disc 
is supported by the high observed luminosity (~ 5 x 10^* 
erg s~^), unachievable for stellar wind accretion. Further- 
more, the pulse period changes at a high rate during the 
outburst (Tsygankov et al., in preparation), which is typi- 
cal for binary systems with disc accretion and observed in 



other transient X-ray pulsars with Be companions (see e.g. 
Tsygankov & Lutovinov 2005 for KS 1947+300). 



6.1 Cyclotron Line 

It was found by iBasko fc SunvaevI ^^76a) that there is a 
critical value of the luminosity (L* ~ 10"^^ erg s~^) dividing 
two accretion regimes: the regime when the influence of the 
radiation on the falling matter is negligible and the regime 
when this influence is significant. When L < L* , the matter 
free-fall zone is extended almost down to the surface of the 
neutron star. In the opposite case (L > L*), observed for 
V0332+53, the radiation-dominated shock rises high above 
the neutron star surface. Almost all of the kinetic energy of 
the infalling gas is lost in this shock, and is then emitted 
late rally by the sides o f the ac cretio n column. 

teasko fc SunvaevI l^Tg^) and iLvubarskii fc SunvaevI 
il98a) showed that the height of the shock H depends on 
the source luminosity 



where m is the dimensionless accretion rate in units 
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Figure 8. Dependence of the ratio of energies of second and main 
harmonics of the cyclotron line on the source luminosity (3-100 
kcV). Triangles and squares represents INTEGRAL and RXTE 
results, respectively. 



of 10^^ erg s~^, 77 - a function depending on the magnetic 
field on the neutron star surface Bms and the thickness of 
the accretion column. As can be seen from the equation, 
the height H changes practically linearly with m in a wide 
range of values, i.e. the shock height grows linearly when the 
source luminosity is increased. 

It was shown above that the energy of the cyclotron 
line detected in the V0332-I-53 spectrum grows approxi- 
mately linearly with decreasing source luminosity. The max- 
imum relative change of the energy and, consequently, the 
corresponding magnetic field are about ~ 25%. In an ap- 
proaching of the dipole field of the neutron star, it corre- 
sponds to a 7.5% relative change of the height h where the 
feature is formed. At the end of the outburst, the source 
luminosity falls to ~ 10^^ erg s~^, the shock descends, 
the column height decreases and we receive emission com- 
ing virtually from the neutron star surface. Because of the 
smallness of the relative change H/Rms, we can consider 
B{h) oc Bms — oih to a first approximation, where a is a co- 
efficient. Comparing with the relation between Ecyci,\ and 
L determined earlier we obtain h oc L. Thus the height h, 
where the cyclotron feature is formed, has the same lumi- 
nosity behaviour as the shock h eight H. 

According to Basko & Sunyaev] l|l976i), only a small 
fraction of the energy accumulated in the accreting matter 
is emitted at the shock. Its main part goes into the extended 
sinking zone below the shock, being gradually emitted by the 
side walls of the accreting column. Moreover, the column 
cross-section has a near-triangular shape where the shock 
proper occupies only a small part in the column central 
plane. In other parts the de celeration of the infalling mat - 
ter is cased by a friction force iLvubarskii fc SunvaevHl98d) . 
Thus, the registered emission is a superposition of emissions 
from different heights above the neutron star surface. There- 
fore, we can consider the height h as some averaged or "ef- 
fective" height of the formation of the cyclotron feature not 
coinciding with the position of the shock itself. 

As seen from Fig|3 the behavior of the energy of the 



second harmonics is qualitatively similar to the main one. 
But due to lower statistics of the data at high energies the 
exact determination of its parameters is model-dependent 
(see above). Thus at the moment we cannot make a final 
conclusion about the rate of change of the second harmonics 
energy with the source luminosity. Nevertheless it is inter- 
esting to note that the ratio of the energies of the second 
and main harmonics slightly decreases with decreasing lu- 
minosity and is approximately equal to the harmonic one 
(2 : 1) near the luminosity ~ 2 x 10^** erg s~^ (Fig|Hl. 

6.2 Pulse Profile 

As experimental information about pulse profiles and their 
evolution has been accumulated, it has become clear that a 
simple model explaining pulsations by the presence of two 
bright spots on the neutron star surface cannot clarify the 
variety of observations. Our results on V0332-I-53 confirm 
this statement. 

As was shown above, the matter flowing from the accre- 
tion disc forms near the magnetic poles accretion columns 
elongated along the magnetic lines of force. Since the falling 
matter is opaque, the radiation is entrained in the col- 
umn and moves with the matter downward, diffusing to 
the edges of the accretion c hannel and escaping laterally 
iLvubarskii fc Sunvaev^^l988^ . Thus a fan beam configura- 
tion o f the X-ray emission is ex pected to prevail in the bright 
state iBasko fc Sunvaevll976al) . which explains the observed 
double-peaked structure of the source pulse profile. But it is 
worth noting that although the radiation has time to escape 
from the column, it will beam toward the neutron star sur - 
face due to relativistic effects llLvubarskii fc Sunvaev^[l988^ 
and simple considerations about the beaming can be inappli- 
cable. Most likely the observed pulse profile variability can 
be explained by a combination of geometrical and physical 
effects. We point out some of them below. 

One of the possibilities to describe the observed changes 
of relative intensities of peaks in the pulse profile can be the 
mechanism of pulse formation f o r large source luminosities 
proposed bv lBasko fc SunvaevI (j[l976b!). Its essence is that 
the magnetic field of the neutron star makes the gas to fiow 
off to the magnetic funnels along the Alfven surface. This 
flow will cover only a part of the surface. The layer of mat- 
ter on the surface will spin with the same angular velocity 
as the neutron star and will periodically shield from the 
observer different parts of emission regions. For a certain 
orientation of the system to the observer it is possible to 
expect different relative intensities of peaks in a dependence 
on the energy band. Also the reflection from the inner sur- 
face of the gas layer flowing to the lower magnetic field can 
make a contribution to t he formation of the pulse profile 
jBasko fc Sunvaev^ll976b^ . When the source intensity is de- 
creased, the scattering optical depth of the masking layer 
is also decreased, which can result in a larger amplitude of 
the asymmetry in the pulse profile and to an appearance 
of some new features in them. The calculations show that 
in bright states the flows on the magnetosphere are opti- 
cally thick. This results in additional changes of the pulse 
profile, especially strong i n the soft and hard energy bands 
llBasko fc Sunvaev|[T976W . 

The most interesting and difficult for explanation 
changes occur near the main harmonics of the cyclotron 
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Figure 9. 3D evolution of the source pulse profiles for 272 (a) and 284 (b) revolutions. Red and blue stripes represents lower and upper 
wings of the cyclotron line. 2D-distributions of pulse profile intensities with the energy and pulse phase are shown in the bottom panel 
by different colors. Solid lines represent levels of equal intensity. Positions of the cyclotron line center are shown by dashed lines. 



line at luminosities lower ^ 7.3 x 10^^ erg s~^. The ob- 
served behavior is difhcult to describe in detail within the 
framework of current models. In partial, it can be con- 
nected with peculiar ities of the radiation beaming near the 
cyclo tron frequ ency llGnedi n fc Sun yaevlll973l : IPavlov et alJ 
ll98fJl . As e.g. iMeszaro^^^geJ I^198!^^ showed, the cy- 
clotron line shape demonstrates a strong angular depen- 
dence. The plasma is more transparent at large angles than 
at small ones for energies below and above the line energy. 
Therefore photons will escape predominantly in the direc- 
tions of large angles, i.e. the radiation beaming in different 
energy channels near the cyclotron line will be strongly dif- 
ferent. In addition, a bulk motion towards the neutron star 
surface also produces a doppler shift, which should result 
in an angular dependence of the cyclotron line energy on 
the viewing angle, that in its turn can give a contamination 
to the observed line energy and pulse profile dependencies 
iBrainerd fc Mezarod Jl99]ji . 



For the better visual perception and understanding 
of the changes described above we built three-dimentional 
pulse profiles with the distribution of their relative inten- 
sities along the pulse phase and energy. To obtain a more 
or less smooth picture we choose the energy window (the 
energy band for each profile) of 4 keV and reconstructed 
a number of pulse profiles with a step of 1 keV from 6 up 
to 45 keV. Such an approach allowed us also to sew to- 
gether JEM-X and ISGRI results around 20 keV relatively 
well. Two such 3D pulse profiles for 272 and 284 revolutions 
are shown in Fig|^ (upper panel). Clearly seen all described 
in section 5 features for both of them. The red and blue 
stripes represents regions of lower and upper wings of the 
cyclotron lines. In bottom panels of Figl^ two-dimentional 
distributions of pulse profile intensities are demonstrated by 
different colors and levels of equal intensities. It is interest- 
ing to trace changes of the maximum intensities for both 
observations: in the first case positions of both peaks are 
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practically unchanged with the energy; in the second one 
the position of the maximum is changed drastically with the 
energy especially near the cyclotron line. The single-peaked 
and double-peaked distribution of intensities in the lower 
and upper wings of the cyclotron line are obviously seen. 

The case of small luminosities is of especial interest for 
following investigations. In this case the radiating plasma 
have an appreciable optical depth only near the cyclotron 
line and the source pulse profile reflects physical properties 
of the plasma flow in the magnetic field, i.e. is determined by 
the anisotropy of the emission and scattering in the plasma. 



7 SUMMARY 

We presented results of the analysis of the INTEGRAL and 
RXTE data obtained during the outburst from the X-ray 
pulsar V0332-I-53. The most important are: 

- for the first time we studied in detail the evolution of 
the cyclotron energy with the source luminosity and showed 
that it is linearly increasing with the source luminosity de- 
creasing in the same way as the change of the height of the 
accretion column; 

- the behavior of the second harmonics energy is qual- 
itatively similar to the main one, but more accurate obser- 
vations are needed for exact measurements of its rate and 
understanding of the behavior of the ratio of energies of sec- 
ond and first harmonics; 

- the strong pulse profile variations with luminosity, 
especially near the cyclotron line, are revealed. 
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